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Abstract 



I discuss a scheme to match perturbative and non-perturbative physics with 
power accuracy in the heavy-quark effective theory. I elaborate on two impor- 
tant aspects of the scheme: 1) a multi-loop subtraction of soft contributions 
from the perturbation series in the pole mass, 2) a perturbative regularization 
of the linearly-divergent heavy-quark self-energy in the lattice formulation of 
the heavy-quark effective theory. 
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In the last few years, the heavy-quark effective theory HQET, has been widely used 
to study hadron physics involving heavy flavors. The HQET is essentially based on the 
concept of scale separation — the large-momentum-scale physics associated with the heavy- 
quark mass can be factorized and calculated in perturbation theory, and the remaining low- 
momentum-scale, non-perturbative physics can be summarized in the heavy-quark effective 
lagrangian. Since the number of papers on the subject is enormous, I refer the reader to 
review articles for relevant references |TJ. 

It was generally accepted that the heavy-quark mass used in the heavy-quark expansion 
should be the pole mass, defined according to the single-particle pole in the perturbative 
heavy-quark propagator 0. Recently, however, Bigi, Shifman, Uraltsev, and Vainshtein 
(BSUV) |§ and Beneke and Brown Q pointed out that the pole mass is intrinsically am- 
biguous at the order of Aqcd and thus, to this degree, the heavy-quark expansion is not 
unique. The discovery is consistent with an earlier observation by Falk, Neubert and Luke 
H that a residual mass term —Sm Q h v h v is generally present in the heavy-quark effective 
lagrangian. It is also consistent with the finding by Maiani, Martinelli and Sachrajda || 
that the lattice formulation of the HQET at order I/tuq is plagued by the divergences in 
inverse powers of the lattice spacing a. 

Given the freedom of constructing the HQET, it is important to find a scheme to match 
perturbative and non-perturbative calculations with power accuracy. Recently, two ap- 
proaches have appeared in the literature to tackle this issue. After arguing that a pertur- 
bative subtraction of power divergences is not viable, Martinelli and Sachrajda JFj proposed 
to use non-perturbative renormalization conditions to define power-divergent operators [§. 
They provided examples to show how it works in the lattice simulations of the HQET. In 
their approach, the cancellation of the infrared renormalons in perturbation series involves 
non-perturbative physics. On the other hand, BSUV [^J took Novikov et al.'s approach to 
Wilson's operator-product expansion (OPE) |§. According to that, the OPE provides a sep- 
aration of physics above and below a certain momentum scale \x. The low-momentum physics 
is taken into account by the non-perturbative matrix elements and the high-momentum 
physics is included in the coefficient functions. In their paper, BSUV showed explicitly how 
this can be done at one-loop level. 

In this paper I follow the proposal by BSUV and aim to complete the scheme in two 
important aspects: I). Outline a procedure for the multi-loop subtraction of the soft contri- 
butions present in the perturbation series in the pole mass. The subtraction can be calculated 
perturbatively using the heavy-quark effective lagrangian in the dimensional regularization 
and minimal subtraction scheme (MS). The approach works also for other perturbation 
series in the HQET. 2). Show how the non-perturbative corrections may be computed, con- 
sistently with the above perturbative calculation, in the lattice formulation of the HQET. 
This involves matching the non-perturbative quantities in the MS and lattice schemes and 
formulating the renormalization conditions for power-divergence quantities in perturbative 
theory. The discussion follows closely a recent paper on Wilson's expansion with power 
accuracy by the present author |jT0|] . 

I start with the heavy-quark expansion for the projected inverse heavy-quark propagator 
in the full theory, as studied by Beneke and Braun |[J, 
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Sp^fc, m Q ) =m Q - m pole (m(/x), a a (fi)) + C( , ct 8 (fi))S~£(v -k,n) + 0{ 5 — , -^-) , 

P m Q m Q 

(1) 

where ttlq, to be specified shortly, is the heavy-quark mass defining the expansion, v is the 
velocity of the heavy quark, and k is the residual momentum. Their relation with the total 
momentum of the heavy quark, P, is 

P = m Q v + k . (2) 

The pole mass m po i e is a function of the renormalized MS mass m(/i). It is also a perturbation 
series in a s , plagued by the infrared renormalon at b = 2ir/f3 in the Borel plane where 
Po = 11 — 2Nf/3. S~ff is a propagator defined in the effective theory with the effective 
lagrangian 

C eS = h v iv ■ Dh v + £ light , (3) 

in the MS scheme, where h v is the heavy-quark effective field and Aight is the lagrangian for 
light quarks and gluons. is a non-perturbative quantity with the ultraviolet renormalon 
at b = 2tt//3q in the Borel plane ||. 

It was suggested in Refs. MM that to cancel the infrared renormalon, one may subtract 
a residual mass 5m(A) of order Aq CD from the pole mass m po i e and choose, 

m Q (A) = m po ie - Sm(A) + ... , (4) 

as the expansion parameter, where the ellipsis represents higher-order terms in the heavy- 
quark mass. Apart from the role of cancelling the renormalon, 5m(A) is otherwise arbitrary. 
In Ref. 0, Sm(A) is identified as the Coulomb energy of the heavy quark below the mo- 
mentum scale A. If so, the expansion parameter mq(A) contains only the physics above 
the scale A and is entirely perturbative if A ^> Aqcd- The A dependence in mg(A) shows 
explicitly that the construction of the heavy-quark expansion is not unique. The scale A 
needs not to be the same as the renormalization scale /x introduced in the full theory. 

To determine the subtraction to all orders, I consider the heavy-quark propagator S~g(v- 
k) in perturbation theory. [I will neglect the generic renormalization scale /i below.] In 
dimensional regularization, the propagator satisfies the renormalization condition, 

S- s \vk = 0, MS)| pert = 0, (5) 

i.e., the mass shift vanishes. The argument for this is simple: The effective lagrangian C c g 
contains no mass scale in the chiral limit. Because 5*^ (v ■ k — 0)| per t has the dimension 
of a mass, it must diverge linearly. All power divergent integrals are taken to be zero in 
dimensional regularization. 

However, the conclusion is deceptive. The self-energy of the heavy quark contains both 
soft and hard contributions, and they cannot cancel each other on the physical ground. 
For definiteness, let us define the hard and soft contributions as follows. Denote a generic 
light-quark or gluon propagator by D(k). Split the propagator into two parts according to 
virtuality k 2 of the particle, 
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D(k) = D{k)6{k 2 > A 2 ) + D(k)9(k 2 < A 2 ) 



(6) 



Then the contribution to the heavy-quark self-energy from a Feynman diagram with all 
propagators replaced by D{k)0{k 2 > A 2 ) is defined as hard. The complementary part, 
denoted as AS, is soft, 

AS(A) = S- S \v ■ k = 0)| pert - S£{v ■ k = 0)| pert (aU k 2 > A 2 ) . (7) 

It is easy to see that AE does not have any linear divergences. All logarithmic divergences 
in it can be regulated dimensionally and subtracted minimally. Therefore we define the 
residual mass, 

5m(A) = C(a s )AS(A) , (8) 

which is a perturbation series in a s , A^ n=1 b n a n s . The one-loop calculation in Ref. || gives 
Ci = 2/3. Multi-loop subtraction can now be routinely calculated according to Eqs. (0) 
and (^). Since 5m(A) contains all the soft contribution, it has the infrared renormalon at 

b = 27T/A,. 

Following Refs. ||(§], I reorganize Eq. (1) as follows, 

Sp l (k,m Q ) = m Q - [m polc - 6m(A)} + C(a s ) [S^(v ■ k, MS) - AS(A)] + ... . (9) 

The infrared renormalon in 6m(A) cancels both the infrared renormalon in m po i e and the 
ultraviolet renormalon in S e s, rendering m po i e — 5m(A) and S~$ — AE(A) well-defined. 
However, both quantities now depend on the separation scale A. 

I now turn to the non-perturbative part of the expansion in the lattice formulation 
of QCD, because at present the lattice provides the only formalism to calculate non- 
perturbative physics from the fundamental lagrangian. In a lattice, the discretized version 
of £ e f[ produces a linearly-divergent quark self-energy. To calculate in the MS-scheme, 
we first match it with the lattice propagator, 

S~g(v ■ k, MS) = Z(fia, a s ) S~$(v ■ k, latt) - ^(O, latt) pert , (10) 

where the second term, calculated in the lattice perturbation theory, diverges like 1/a and 
serves to cancel the linear divergence in the first term. Z is a renormalization constant, 
calculable in perturbation theory fTT||. Alternatively, one can add a residual mass term 



—Smoh v h v to the original lagrangian in Eq. (|3|), where 8mo is fixed by the perturbative 
renormalization condition, 

S£(0, latt)| pcrt . = . (11) 
The propagator calculated in the new effective theory has no linear divergences. The new 



propagator replaces the two terms in the right-hand side of Eq. (10). However, the residual 
mass thus determined, 5m = l/a^ n c n a", has the infrared renormalon at b = 2n//3 . 

There is yet another way to define the effective theory on the lattice, which is free of both 
the linear divergence and the infrared renormalon. Define a new SrriQ with a subtraction of 
another series, 
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5m = 1/a J2 c na n s ~ A J2 d ^ , ( 12 ) 

n=l n=l 

where d n is adjusted so that the following renormalization condition is satisfied, 

S£(0, latt)^ = -Z~ l AZ(A) . (13) 

Then 8m is free of the infrared renormalon, although it now depends on the separation 
scale A. The modified effective lagrangian, 

£'eff(A) = £ cff - 5m (a, A)h v h v , (14) 

is also A dependent, as discussed in Refs. [f§|jj. The non-perturbative heavy-quark propa- 
gators S~fi(v ■ k,A) in this theory, apart from the renormalization factor Z, gives directly 
S^f(v ■ k, MS) — AE(A) in Eq. (EH). Let me emphasize again that the residual 5m is an 
entirely perturbative quantity. 

Using the expansion parameter defined in Eq. (§), I have, 

S P 1 (k,m Q )=C(a s )ZS- s 1 (vk,A) + ... . (15) 

It is possible to find a special value of A = A such that, 

5S(0,Ao)=0. (16) 

The expansion parameter ttiq define at this Ao coincides with the subtracted pole mass 
introduced by Martinelli and Sachrajda. In some sense, this special mass is the physical 
pole mass because the propagator in the full theory has a pole at j>= tuq, as is clear from 
Eq. fllSD . ttIq is independent of the details of the effective theory. It would be nice to check 
that the two formulations of the lattice renormalization conditions give the same w,q. 

Both the subtracted pole mass and the physical mass of a pseudo-scalar meson with the 
heavy quark can be expanded in the mg(A), 

m s Q = m Q {A) + C(a s )S; s \0, A) + ... , 

m H = m Q (A) + C'(a s )A(A) + ... . (17) 

A is the mass of the pseudo-scalar meson calculated in the lattice effective theory with the 
renormalization condition in Eq. (|13|) . C and C are the coefficient functions calculable 
in perturbation theory. A dependence is explicit in different orders of tuq in the above 
equations, however, it cancels in the physical quantities. 

In conclusion, I have discussed in this paper the heavy-quark effective theory in the ap- 
proach of BSUV. I gave a multi-loop subtraction formula for calculating the scale-dependent 
heavy-quark mass 7tiq(A). I outlined a perturbative regularization of the linear divergences 
in the heavy-quark self-energy in the lattice QCD. With these crucial ingredients, power 
corrections to the results of the leading-order HQET can be practically calculated. 
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